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ABSTRACT 

We report on XMM-Newton spectroscopic observations of the luminous, radio- 
quiet quasar PDS 456. The hard X-ray spectrum of PDS 456 shows a deep absorption 
trough (constituting 50% of the continuum) at energies above 7 keV in the quasar rest 
frame, which can be attributed to a series of blue-shifted K-shell absorption edges due 
to highly ionized iron. The higher resolution soft X-ray grating RGS spectrum exhibits 
a broad absorption line feature near 1 keV, which can be modeled by a blend of L-shell 
transitions from highly ionized iron (Fe XVII - xxiv). An extreme outflow velocity 
of ^ 50000 km s‘^ is required to model the K and L shell iron absorption present 
in the XMM-Newton data. Overall, a large column density (N^ = 5 x 10~^cm~~) of 
highly ionized gas (log ^ = 2.5) is required in PDS 456. A large mass outflow rate of 
^ 10 M© year"* (assuming a conservative outflow covering factor of 0.1 steradian) is 
derived, which is of the same order as the overall mass accretion rate in PDS 456. This 
represents a substantial fraction (~ 10%) of the quasar energy budget, whilst the large 
column and outflow velocity place PDS 456 towards the extreme end of the broad 
absorption line quasar population. 

Subject headings: galaxies: active — quasars: individual:PDS456 — X-rays: galaxies 


1. Introduction 

PDS 456 is a luminous, but low redshift (z = 0.184) radio-quiet quasar identified in 1997 
(Torres et al. 1997). The optical and infra-red spectra (Simpson et al. 1999) show broad Balmer 
and Paschen lines (e.g. H/3 FWHM 3000 km s"*), strong Fe ll, a hard (de-reddened) optical 
continuum (/i, oc *^^ *), and one of the strongest ‘big blue bumps’ of any AGN (Simpson et al. 
1999, Reeves et al. 2000). PDS 456 has a de-reddened, absolute blue magnitude of M^ w ~27 
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(Simpson et al. 1999), making it at least as luminous as the radio-loud quasar 3C 273. Indeed PDS 
456 is the most luminous known AGN in the local Universe (z < 0.3), its luminosity being more 
typical of quasars at z=2-3, at the peak of the quasar luminosity function. 

PDS 456 was first detected as the X-ray source RXS J172819.3-141600 in the ROSAT All Sky 
Survey (Voges et al. 1999). Subsequent ASCA and RXTE and BeppoSAX observations showed that 
it was highly X-ray variable (Reeves et al. 2000, 2002), requiring that the X-ray source was very 
compact (< 3 Schwarzschild radii) and accreting with a unusually high efficiency (between 6% 
and 41%), close to the limit for a maximal Kerr black hole (Thome 1974). The X-ray spectmm 
of PDS 456, obtained by ASCA and RXTE also appeared to be unusual for a quasar. A strong 
edge like feature was present in the iron K-shell band, either due to a highly ionized absorber or 
reflector, whilst absorption towards the soft X-ray band was also observed (Reeves et al. 2000). 
An XMM -Newton AO-1 observation of PDS 456 was subsequently scheduled, in order to obtain 
high signal to noise spectra of both the iron K-shell complex with the EPIC detectors and a high 
resolution spectrum of the soft X-ray absorber with the XMM-Newton RGS (Reflection Grating 
Spectrometer). We report on the X-ray spectra obtained through the XMM-Newton observations 
of PDS 456, which reveal evidence for a massive outflow, with a velocity of ~ 50000 km s"‘. 
A previous paper (Reeves et al. 2002) discussed in detail the X-ray variability during the XMM- 
Newton observation and a simultaneous week long BeppoSAX observation. 


2. The X-ray Observations 

PDS 456 was observed by XMM-Newton on 2001 February 26th with a duration of 40 ks. 
Data were taken with the EPIC-pn detector (Struder et al. 2001) in Full Window Mode and with 
the EPIC-MOS detector (Turner et al. 2001) in Large Window mode. The data was reduced using 
version 5.3 of the XMM-SAS. All time intervals were included, as the observation contained no 
background flares. Data were selected using event patterns 0-12 (for the MOS) and pattern 0-4 
(for the pn) and only good X-ray events were included. The spectra were extracted from circular 
source regions of 40" radius, whilst background spectra were extracted from an offset circle, close 
to PDS 456, but free of any background sources. Response matrices and ancillary response files 
were generated using the SAS tasks RMFGEN and ARFGEN respectively. RGS spectra were also 
reduced using the pipeline SAS script RGSPROC. 

Background subtracted spectra were fitted using XSPEC vl 1.2, including data over the energy 
ranges 0.4 to 12 keV for the EPIC-pn and 0.4 keV to 10 keV for the co-added EPIC-MOS spectrum. 
A Galactic absorption column of = 2 x 10^' cm“^ (Dickey & Lockman 1990) was included in 
all the fits. A single absorbed power-law fit (photon index, F = 2.0) to both MOS and PN data 
is clearly inadequate (reduced = 10). Figure 1 shows the data/model ratio to this power-law, a 


-3- 


strong soft excess is observed below 1 keV and a large drop in counts is observed above 7 keV 
in the iron K-shell band. The MOS and pn spectra are consistent in spectral shape, except below 
0.5 keV where there is a small divergence in spectral slope. 


3. The Ionized Absorber in PDS 456 

Initially, to model the iron K-shell band, the data were fitted from 2-12 keV (or 2-10 keV 
for EPIC-MOS), so that the effects of the large soft excess are negligible. As the EPIC-pn data 
contained significantly more counts than the EPIC-MOS at the highest energies, we proceeded to fit 
the EPIC-pn data only. The ratio spectrum above 2 keV to a power-law fitted over the 2-5 keV range 
(i.e. outside of the iron K band) is shown in the inset to Figure 1. A strong deficit of counts occurs 
above 7 keV, where one would expect to observe K-shell absorption edges from highly ionized 
iron. Indeed a simple Galactic absorbed power-law produced a poor fit (x'/^of = 1001/746, null 
hypothesis probability 9 x 10”*®). 

Adding three absorption edges between 7 and 10 keV significantly improves the fit ix^fdof= 
739/740), modeling the broad Fe absorption feature present in the X-ray spectrum. The energies 
(in the quasar rest frame) and optical depths of the edges are then; £ = 9.3 ± 0. 1 keV, £ = 8.5 ± 
0.2 keV and £ = 7.6±0.3 keV withr = 0.56±0.11, r = 0.26±0.11,t = 0.24±0.08 respectively. 
The breadth of the absorption feature suggests that a wide range of ionizations may be present, 
whilst the depth of the feature implies a column of gas of Nh ~ 10*^ cm”". Note that the upper 
limit to a narrow (i.e. unresolved) Fe emission line (in the energy range 6.4 - 6.97 keV) is < 12 eV, 
implying that the solid angle subtended by the absorber must be small (Q < 0.5 steradian). Note 
that an extreme broad iron Ka emission line from an accretion disk cannot model the shape of 
the spectrum from 2-12 keV. If one attempts to model the spectmm with a Fe line profile from an 
extremal Kerr black hole (e.g. Laor 1991), the physical parameters required are unphysical (disk 
emissivity inclination=90“, equivalent width > 1 keV), because the line is (unsuccessfully) 
attempting to model the broad absorption dip above 7 keV with the blue-wing of the line profile. 

To obtain a more physical representation of the Fe K-shell absorber, we used a grid of models 
generated by the XSTAR photoionization code (Kallman et al. 1996) to fit the EPIC-pn data above 
2 keV. Solar elemental abundances were assumed, whilst a turbulence velocity of 1(K)0 km s”* was 
used. A one zone photoionization model was adopted, with the column density (^Vh), ionization 
parameter (^), partial covering fraction (/) and outflow velocity (vout) as free parameters. An 
effective hydrogen column density (for a solar abundance H/Fe ratio) of 5.71^5 ^ ‘ is 

required, together with an ionization parameter of log ^ = 2.5 ± 0.3 and a partial covering fraction 
of ~60%, whilst the power-law continuum steepens to F = 2.40 ± 0.05. At this ionization state 
most of the absorption is due to Fe xvii - Fe xxiv (i.e. where the L-shell is partially filled). The 
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best fit spectrum modeled by XSTAR is shown in Figure 2 (panel a), the depth and energy of the 
absorption feature being well matched, whilst the partial covering warm absorber also models the 
convex spectral curvature (through the increase in opacity at lower energies due to iron L-shell 
absorption) present in the EPIC-pn data between 2-6 keV. Overall, the fit statistic for this model is 
excellent (x~ldof = 717/740). 

Importantly a large outflow velocity of 47000 km s”' (90% confidence limits, 36000 km s"* 
- 82(KK) km s"') is needed to model the spectrum. Although the exact outflow velocity is not well 
constrained by the EPIC-pn data alone, a significant outflow velocity is still required. Figure 2(b) 
illustrates this point, if one fixes the absorber velocity to zero, then the model fit becomes poor as 
the iron K edges then lie at a significantly lower energy than is observed in the EPIC data. Indeed if 
one tries to increase the energy of the model absorption feature by increasing the ionization state of 
the absorber (so that He and H-like iron dominates), the result is a worse fit still, partly because the 
absorption edge feature becomes weaker (as a significant fraction of iron becomes fully ionized) 
and also because the most prominent absorption features are now the strong Ka and K/5 absorption 
lines due to Fe XXVl which are not present in the actual XMM-Newton data (Figure 2c). Thus a 
significant outflow velocity is required by the data, indeed one can place a conservative lower-limit 
(at 99.9% confidence) on the velocity of 24000 km s”' (or ~ 0.08c). 


3.1. The soft X-ray RGS Spectrum of PDS 456 

Extrapolation of the hard X-ray spectrum to lower energies (down to 0.4 keV), leaves clear 
residuals in the EPIC MOS and pn data (e.g. Figure 1). A strong soft excess is present below 1 
keV, whilst residuals due to an additional absorber may also be present in the soft X-ray band. To 
investigate this at higher spectral resolution, we used data available from the Reflection Grating 
Spectrometer (RGS) below 2 keV. The most noticeable feature present is a broad absorption trough 
observed between 12-15 A (see Figure 3). The feature is highly significant, being detected at 
> 99.99% confidence. When parameterized by a simple broad Gaussian shaped absorption line, 
the best fit parameters (in the quasar rest frame) are E = 1.07 ± 0.02 keV with an equivalent width 
of 26 ± 10 eV and a velocity width of 9000 d= 3000 km s"' . 

As in the previous section, we used our grid of photoionization models to characterize the 
properties of the absorbing gas, apparent in the RGS spectrum. A lower colunm density (iVn ~ 5 x 
10" cm'-) is required to fit the data, whilst both a similar ionization parameter (log ^ = 2.6 ±0.5) 
and outflow velocity (with Vout = 570001*o^ km s'*) are needed, as per the fits to the EPIC-pn 
data. Partial covering is not required in the RGS fits. The photon index obtained (F = 2.4, after 
correcting for absorption) is consistent with the EPIC-pn data, whilst a soft blackbody component 
(with kTsB = 80 ± 6 eV) is also needed in the fit. The lower column density derived from the RGS 
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fit indicates that the absorber is patchy and thus is likely to be located quite close to the X-ray 
source. Thus we see a higher (partially covered) column along the line of sight towards the hard 
X-ray continuum, whilst the column towards the soft X-ray continuum is lower. 

Inspection of the best-fit photoionization model implies that a blend of L-shell transitions 
from Fe XVII - Fe XXIV as well as K-shell absorption fi'om Ne ix/x is responsible for the broad 
absorption feature at 12-15 A. A weak Ly-a absorption line due to O vill may also be present at 
19 A. The outflow velocity derived appears robust, due to the positioning of the broad absorption 
feature. If one fixes the absorber velocity at zero in the quasar rest-frame, then the Fe L-shell ab- 
sorption in the model is located at energies below the observed absorption trough, whilst increasing 
the ionization of the absorber only serves to significantly weaken the soft X-ray absorber (as all 
the iron L-shell electrons are then removed from their ions). 


4. Discussion 

X-ray spectroscopy of PDS 456 with XMM-Newton has revealed a high column density (Ah = 
5 X 10'^ cm~-), highly ionized (log ^ = 2.5) warm absorber, its main signature being a series of 
K and L shell absorption edges and lines due to iron in the ionization range Fe XVll - Fe xxiv. 
Crucially, a near relativistic outflow velocity is required, of the order 50000 km s~’ . This velocity 
is similar to the actual systemic velocity of PDS 456 (at z = 0.184), which means that it is possible 
that the absorber is actually located nearby at z = 0 (i.e. with no net outflow). However we consider 
it extremely unlikely that such a high column, highly ionized absorber could occur by chance along 
the local line of sight. Furthermore, although the Galactic column towards PDS 456 is moderately 
high (Ah = 2 x 10"'cm"'), this neutral absorption does not contribute towards the features seen 
in the iron K and L-shell bands, and in any event has been accounted for in all our modeling. 
Therefore we consider it highly probable that the absorber is local to PDS 456 itself. 


4.1. A Massive Outflow in PDS 456 

The combination of a large column density, high ionization and a fast outflow velocity in 
PDS 456 indicates that the absorber is located close to the quasar central engine and that its mass 
outflow rate is unusually large. A maximum distance to the absorber can be calculated on the 
condition that its thickness (A/?) cannot exceed its distance (/?), i.e. A/?//? < 1. This requires that 
R < 10^/?g (where Rg = GMlc~) for a black hole mass of IO^A/q for PDS 456, thus it is located 
within the BLR. One likely possibility is that the absorber arises through a disk driven outflow or 
wind. The exact origin of these winds and the mechanism responsible for driving them (radiation 
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and/or magnetic) is unclear. For Broad Absorption Line (BAL) QSOs, line driven winds have been 
proposed (e.g. Murray et al. 1995, Proga, Stone & Kallman 2000) as well as centrifugal driving 
along open magnetic field lines (e.g. Bottorff et al. 1997). 

Using conservation of mass, the outflow rate (Mow) for an outflow subtending a solid angle of 
n steradian with a constant velocity Vout is given by> 

Mow = finr^VoutWp 

where n is the density, r the distance from the X-ray source and nip the proton mass. Now the 
ionization parameter ^ is defined as ^ = Lx/nn, so we can obtain:- 

Mow “ f7vout^pLx/ ^ 

From the observations, the hard X-ray luminosity of PDS 456 (above 2 keV) is ~ 10'*^ erg s'*, 
the ionization parameter is ^ ~ 1000, whilst v = 50000 km s'*. From the lack of any iron K line 
emission it is likely that the outflow subtends a solid angle of < 0.5 steradian. However even 
for a conservative value of 0.1 steradian, the mass outflow rate derived is M<,« = 8 x 10*^ g s'* or 
~ IOMq year"’. An absolute lower limit of ~ 5M© year'* is obtained using the 99.9% confidence 
lower bound on the outflow velocity of Vout = 25000 km s'* . PDS 456 is an extremely luminous 
quasar with a bolometric luminosity of Lbot ~ 10^’ erg s'* and a black hole mass of lO^M© assum- 
ing accretion occurs near the Eddington limit (Simpson et al. 1999; Reeves et al. 2000). The mass 
accretion rate for PDS 456, given by M = Lboi/rjc-, is ~ 20M© year'*, for an efficiency of ?; = 0.1. 
Thus it seems likely the outflow rate is of the same order as the mass accretion rate in PDS 456. 
Similarly the kinetic output of the outflow in PDS 456 is large, about lO'*^ erg s'* corresponding to 
10% of the bolometric output of the quasar. 

The high velocity also implies that the outflow is likely to be driven from the accretion disk 
at a fairly small radius. For a Keplerian disk, the launch radius for the outflow is simply Rgut = 
c-/v\j(g, where Vesc is the escape velocity of the material. Assuming that Vesc ~ Vout (although 
some acceleration of the matter will probably occur), this implies that the outflow in PDS 456 is 
launched from a radius of ~ 40/?g (about 10*^ cm for a lO^M© black hole). 

PDS 456 appears to be in a small but growing group of AGN to show extreme absorption 
properties in the iron K-shell band. The lensed BAL quasar APM 08279+5255 has been found to 
exhibit iron K absorption lines and/or edges (Chartas et al. 2002, Hasinger, Schartel & Komossa 
2002) in both Chandra and XMM-Newton observations. A similar absorber has also been reported 
in another BAL quasar PG 1 1 15+080 (Chartas, Brandt, & Gallagher 2003). The outflow velocities 
in these cases are similar to PDS 456, in the range 0.1c -0.4c. Given this apparent association 
between high velocity X-ray outflowing gas and the UV BALs, one may predict to see BAL like 
features in the UV spectrum of PDS 456. Indeed a HST-STIS spectmm of PDS 456 (O’Brien et 
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al. 2003, in preparation) reveals that there may be cooler gas present in the outflow; a BAL is seen 
in Lyman-a with an outflow velocity of 12000 - 22000 km s"', whilst the C IV emission line is 
also blue-shifted by ~ 5000 km s"'. Thus PDS 456 may just reside at the extreme end of the BAL 
quasar population (e.g. Tumshek et al. 1988; Weymann et al. 1991). 

Similar highly ionized X-ray outflows have recently been discovered in the non-BAL radio- 
quiet quasars PG 1211+143 (Pounds et al. 2003a) and PG 0844+349 (Pounds et al. 20O3b), suggest- 
ing that these properties are not unique (but are probably related) to UV identified BAL quasars. 
We also note that strong iron K absorption has also been observed in the NLSls IH 0707-495 and 
IRAS 13324-3809 (Boiler et al. 2002, 2003) and also in some X-ray binaries (e.g. GRS 1915+105, 
Lee et al. 2002; GX 13+1, Sidoli et al. 2002; X 1624-490, Parmar et al. 2002). In the quasars PDS 
456, PG 121 1+143 and PG 0844+349, the X-ray columns are near to 10-‘*cm”*, the outflow veloc- 
ities measured are ~ 0.1c, whilst the ionizations are also similar (^ ~ 10^). Using PG 1211+143 
as an example. King & Pounds (2003) have recently demonstrated, for radiatively driven out- 
flows, that objects accreting near to the Eddington rate are likely to contain massive (Afom ~ Afedd) 
Compton-thick winds. Indeed when the outflow rate is of the same order as the accretion rate (as 
is the case in PDS 456), the flow can become optically thick at ~ 10- 100/?j. This can account for 
both the high columns and velocities of the X-ray outflows, as well as provide a natural source of 
thermalized photons for the UV and soft X-ray excesses observed in these quasars. 

However other scenarios are possible, for instance Reeves et al. (2002) recently used a mag- 
netic flare model to explain the rapid X-ray variability in PDS 456, through coherent flaring events 
occuring when the accretion rate is close to Eddington. Indeed it is conceivable that the outflowing 
matter in PDS 456 is the ejecta associated with these flares. Interesting the output of the magnetic 
events (~ 10'*^ erg s"') in PDS 456 is very similar to the kinetic power of the outflow, so this sce- 
nario seems energetically plausible. Overall the extreme properties of the outflow in PDS 456 seem 
linked to its high accretion rate (e.g. Reeves et al. 2000), which appears to drive the energetics of 
its unusually massive outflow. 
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Fig. 1. — Ratio of the broad band X-ray spectrum of PDS 456 to a simple Galactic absorbed 
power-law model (with F = 2.0). EPIC-pn data are plotted in black, EPIC-MOS in gray. A clear 
soft X-ray excess is observed below 1 keV, as well as a large decrease in counts above 7 keV. The 
inset plots the data/model ratio to a power-law fit to the 2-12 keV EPIC-pn spectrum in greater 
detail. A 50% drop in counts above the iron K-shell band is observed. Note the energy scale is 
plotted in keV in the observed frame. 
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Fig. 2. — Warm absorber fits to the EPIC-pn data of PDS 456 using the XSTAR photoionization 
code, (a) shows the best fit model (solid line) fitted to the EPIC-pn data (crosses), which requires 
an outflow velocity of ~ SfXKX) km s"’ and an ionization parameter of log ^ = 2.5; (b) shows the 
same model, but with no outflow velocity, resulting in a significantly worse fit; (c) the same model, 
but with no outflow and a higher ionization parameter of log ^ = 3.5. At this high ionization the 
model cannot reproduce the shape of the iron K-shell absorption seen in ft)S 456. 
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Fig. 3. — XMM-Newton RGS-2 spectrum of PDS 456, the data points are shown as crosses. A 
strong absorption trough is present at 12-15 A, observed frame. The best fit photoionization model 
(solid line) is superimposed on top of the data points, a high outflow velocity of ~ 50000 km s"* is 
required for the absorber. Note due to a malfunctioning chip, RGS-1 does not containing data in 
the 12-15 A range and thus is not plotted. Most of the opacity of the absorber results from L-shell 
transistions of iron in the ionization range Fe xvii - Fe XXIV. 



